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The study of QCD at finite baryon density is severely hampered by the so-called fermion sign problem. As a 
result, we have no known first principles approach to study nuclear matter, or neutron stars from QCD. On the 
surface, the large N c limit does not seem to simplify matters. In this limit, however, one can exploit dualities 
that exist between strongly coupled gauge theories. Our focus will be on some rather novel orbifold equivalences 
that have recently been discovered at finite density. These equivalences relate strongly coupled theories plagued 
by a sign problem, to strongly coupled theories free of sign problems. As a result, such dualities give deeper 
insight into the nature of the sign problem and possibly provide a way to simulate QCD at finite density in the 
large N c limit. 



D 
r- 1 

i — i 1 . Motivation and Scope 

Lattice gauge theory provides a quantitative, numerical technique for answering non-pcrturbative questions 
about strongly interacting theories from first principles. As simulations of QCD and other gauge theories are 
carried out in Euclidean space by computational necessity, there is a fundamental restriction on the types of 
non-perturbative questions that can be answered using lattice techniques. In particular, the study of QCD at 
finite baryon density is not possible using Euclidean lattice simulations. This is due to the fermion sign problem. 
In the statistical evaluation of the QCD path integral, one uses the weight 

V[A^] =Det(X>L4 M ])e- s ™^, (1) 

to generate an ensemble of gauge configurations, A^(x), in order to sample the important contributions to the 
• i-h QCD partition function. Notice we employ a continuum formulation throughout for simplicity. The determinant 
factor arises from integrating out fermionic fluctuations, and depends on the fermion action, T>, which has the 
form 

03 

V[A^ = (d M + iA^) ltl + ^74 + to. (2) 

Here we have included a chemical potential coupled to fermion number, /i. 1 At vanishing chemical potential, 
/j = 0, the eigenvalues of come in complex conjugate pairs. The determinant is consequently real, and one 

can use 'PfAJ as a probabilistic weight to generate gauge configurations. On the other hand, when /j, ^ 0, the 
fermion sign problem manifests itself as a complex-action problem: "P[j4 m ] is complex- valued, and thus cannot 
be utilized as a probabilistic weight. 

As a consequence of the fermion sign problem, we have no known first principles technique to address the 
physics of QCD at finite baryon density. For systems with baryon number on the order of a few, one can build 
these light nuclei three-quarks-at-a-time using lattice QCD correlation functions. While there is a signal-to-noise 
problem infecting such computations, impressive progress is being made. Current simulations are now reaching 



baryon numbers of Nb = 3 and 4 for the first time, see Beane et al. 2009 , Yamazaki et al. 2010 . Beyond light 



*Electronic address: a.cherman@damtp.cam.ac.uk 
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1 In QCD, the quark chemical potential, fi q , is trivially related to the baryon chemical potential, /ig, by the relation fig = 3/iq. 
We will often be sloppy and simply refer to the quark chemical potential as the baryon chemical potential. One should keep in 
mind that in SU(N C ) gauge theory, there is a factor of jV c relating the two. 
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nuclei, one wishes to understand how QCD degrees of freedom manifest themselves in larger nuclei (such as 
observed in deep inelastic scattering off nuclei, the so-called EMC effect), and in nuclear matter (such as in the 
passage of jets through dense media in heavy ion collisions). Finally, the study of dense quark matter is likely 
relevant to describe the interior of compact stars. The density in the interior of a neutron star, for example, 
might be great enough to support Cooper pairing of quarks leading to color superconductivity, see | Afford et al. 



2008 , Rajagopal and Wilczek 2000 . Determining whether this is actually the case remains an outstanding 
theoretical challenge. While models used to study high density QCD are becoming more realistic, a definitive 
answer directly from QCD remains elusive. 

Unfortunately, we shall not attempt to solve the QCD sign problem here. Instead, we shall garner insight 
into the problem by studying certain QCD-like theories. A QCD-like theory we take to be any gauge theory 
with fermion degrees of freedom, which does not suffer a sign problem at finite density. For example, it is 
well known that SU{2) gauge theory does not suffer from the fermion sign problem at finite baryon density, 
and has been the subject of numerous theoretical and lattice studies; see, for example, two early works on the 
subject, Hands et al. [f999 , Kogut et al. 1999 . There are various other gauge theories which do not suffer 
from a sign problem at finite density. For gauge groups of arbitrary rank, there is also: SU (N c ) gauge theory 
at finite isospin density, SU(N C ) gauge theory with adjoint matter at finite baryon density, SO(N c ) gauge 
theory at finite baryon density, and Sp(2N c ) gauge theory at finite baryon density (the latter two can also be 
considered at finite isospin density). While these theories are each amenable to Euclidean space Monte-Carlo, 
a natural question emerges: how like QCD are these QCD-like theories? As we shall focus on the large N c limit 
of QCD, we wish to know specifically how these theories are like SU(N C ) gauge theory at finite baryon chemical 
potential. 2 

Surprisingly it can be shown precisely how these theories are related in the large iV c limit, Hana da and| 



Yamamoto] [201 1| . The quantitative relation between these QCD-like theories can be elucidated due to recently 



discovered orbifold equivalences at finite density. The first such equivalence was found in Cherman et al. 2011 



Instead of exploring the universality of QCD-like theories at finite density, we shall focus on just one QCD- 
like theory, namely that of SO(2N c ) gauge theory at finite baryon chemical potential. This theory has yet 
to be studied using lattice Monte-Carlo techniques, but could provide quantitative information about QCD at 
finite density. Our focus throughout is on this connection. We begin by exploring the orbifold equivalence of 
SO(2N c ) gauge theory to SU(N C ) gauge theory at finite baryon density. While the orbifold equivalence holds 
in perturbation theory, it has been shown to break down non-perturbatively as the baryon chemical potential is 
increased. We then investigate whether certain deformations of SO(2N c ) gauge theory can extend the validity 



of the equivalence to larger values of the baryon chemical potential, jCherman et aL |2011 , Cherman and Tiburzi 



2011 . Our presentation ends with a brief summary, and an outlook to future work. 



2. Orbifold Equivalence at Finite Density 



The gauge group SO(2N c ) contains SU(N C ) as a subgroup. Projecting out everything but these SU(N C ) 
degrees of freedom from SO(2N c ) allows us to establish an orbifold equivalence between the two theories valid 
in the large N c limit. As we wish to include matter fields as well, we will spell out the case of a more general 
orbifold equivalence between two theories. The notion of orbifold equivalence was first discovered in the context 
of string theory, Kachru and Silverstein |1998 . Shortly after it was realized t hat such equivalences also applied 



to supersymmetric gauge theories, as well as non-supersymmetric theories, see Bershadsky and Johansen 1998 



Schmaltz 1999 



For our purposes, we spell out an orbifold equivalence in three steps, and specialize to the particular case of 
interest below. 



1. Consider a gauge theory, which we shall refer to as the parent theory 
under a discrete subgroup Zr- 



2. Form what is referred to as the daughter gauge theory 
parent theory P. 



. Suppose that P is invariant 
which consists of all Zr-invariant fields of the 



2 The phenomenological relevance of learning about the behavior of large N c QCD depends on the closeness of the large N c 
world to ours, which is a subtle question discussed, for e xample, in|Buchoff et al.| [2010| , |Frandsen et al.] [2006] , |McLerran and| 



Pisarski 



2007 



Park et al. 



2000 



|Shuster and Son 



2000 



Torrieri and Mishustin 



2010 
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3. There is a planar equivalence between correlation functions of Zp-neutral fields of P and the corresponding 
correlators in D. 

It is useful to spell this equivalence out symbolically. Writing the fields of the parent gauge theory P as <fi® r , 
where Qr is the charge of the field <fip under Zp, and the fields of the daughter gauge theory D as (pa, we can 
write the orbifold equivalence of correlators generically in the form 

(■••tf r=0 ---> Ar ' = O °(---^--->- (3) 

The equivalence generally relates only a subset of the correlation functions that exist in D. 

The idea of using a orbifold equivalence at finite density is to embed the sign problem in the equivalence. 
Thus one hopes to have an equivalence between a theory without a sign problem, and one with a sign problem. 
In the non-perturbative regime, both theories in the equivalence are strongly coupled; however, the parent is 
amenable to lattice Monte-Carlo, while the daughter is not. Hence, one can use the equivalence to compute 
certain correlators of the daughter theory indirectly from the parent. 

Now consider SO(2N c ) gauge theory coupled to matter in the fundamental representation. At any value of 
the baryon chemical potential, the theory has the following discrete symmetry 

Z 2 : 1 ^^e"/ 2 eSO(2N c )<S)U(l) B . (4) 

The daughter theory formed from Z2-invariant operators can be shown to be SU(N C ) gauge theory 3 with 
fundamental matter and a baryon chemical potential, Cherman et al. [2011] . Thus one can compute correlations 



functions of neutral operators in SO(2N c ) gauge theory to learn about large N c QCD at finite baryon density. 
One such correlation function, for example, is the chiral condensate, (f/'V')- A crucial point to note is that 
the discrete symmetry must involve the baryon number transformation in order to project onto SU(N C ) gauge 
theory at finite baryon density. 



3. Equivalence Breakdown 



We have spelled out the orbifold equivalence between a sign-problem-free theory, SO(2N c ) at finite chemical 
potential, and a sign-problematic theory, large N c QCD at finite chemical potential. At this point, it just seems 
we are hiding the fermion sign problem in a complicated theoretical apparatus. In fact, this is not too far from 
the truth. While orbifold equivalences hold to all orders in perturbation theory, there are well-known problems 



with such equivalences at the non-perturbative level, see, Armoni et al. 2005 , Kovtun et al. 2005 , Tong 2003 



Unsal and Yaffe 



2008 



Generally what happens is that phase transitions can spoil orbifold equivalences. For 
example, the parent theory can develop condensation which spontaneously breaks the symmetry involved in the 
orbifold projection. A heuristic picture of why such symmetry breaking is detrimental to orbifold equivalences 



is sketched in Cherman and Tiburzi 2011 



For the specific orbifold equivalence at hand, we must realize that baryon number can be spontaneously 
broken in theories with an orthogonal gauge group. This is because such theories contain two types of baryons, 
the type- A baryons formed of N c fermions, and the type-i? baryons formed of just two fermions, 



The diquark interpolating field, Bb, while gauge invariant, is not invariant under the Z2 symmetry used in the 

orbifold projection, Eq. Q, because B B —Bb- In terms of correlators, we might not be too worried about 
such states since they have no image in the daughter theory. Unfortunately this is naive. Because the Bb are 
charged under baryon number, the chemical potential will lower the energy of such states. Once the chemical 
potential is large enough, it will become energetically favorable for them to condense. The Z2 symmetry will 
then be spontaneously broken. One might further hope that the size of the chemical potential needed for 



s Strictly speaking, the daughter theory actually has the gauge group U(N C ); however, the difference between this gauge group 
and SU(N C ) is suppressed at large N c . 
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Figure 1: Phase diagram of SO(2N c ) gauge theory. The region of spontaneously broken baryon number corresponds 
to the region where orbifold equivalence to large JV C QCD breaks down. Outside this shaded region, the equivalence is 
valid. Other phase transitions (chiral, deconfinement) are not shown. 



such condensation is large; however, it turns out that these states, having scalar diquark quantum numbers, are 
abnormally light. The SO(2N c ) gauge theory with Nf light quarks has an enhanced SU(2Nf) flavor symmetry 4 
that is spon taneously broken down to S O(2Nf) at vanishing chemical potential by the formation of the chiral 



Coleman and Wittcn 



1980 



. The resulting N f (2N f + 1) 



condensate, 

pseudoscalar pions, and Nf(Nf + L) baryonic pions, Cherman et al. 2011 , Cherman and Tiburzi 2011 , Peskin 



1 Goldstone modes consist of — 1 



1980 . The latter turn out to have the quantum numbers of Bb and Bb- The Bb modes will condense at zero 
temperature for a rather small value of the chemical potential, fis = \ m -m and invalidate the equivalence. Thus 
the equivalence is valid for fis < h m-K- 

The phase diagram for SO(2N c ) gauge theory is shown in Fig. [I] where we depict the region of spontaneous 
baryon number violation. For moderate values of the chemical potential, a Bose-Einstein condensate (BEC) of 
baryonic-pions will form. At asymptotically high densities, one can show that Cooper pairs will form in the scalar 
diquark channel. These two regions are very likely continuously connected. 5 Notice that the phase diagram 
for SO(2N c ) gauge theory with baryon chemic al potential is qualitativel y similar to SU(N C ) gauge theory at 
finite isospin density, which was first studied in Son and Stephanov 2001 . In fact, the similarities can be made 



quantitative in the large N c limit, because there is an orbifold equivalence relating the two theories, Hanada 
and Yamamoto 2011 . This orbifold equivalence does not involve a baryon number transformation, and is 



conjectured to be valid in the entire pb-T plane. 



4. Deformations 



We have seen that the orbifold equivalence between SO(2N c ) gauge theory at finite baryon density and large 
N c QCD at finite baryon density only holds in a limited region. Beyond a certain value of the chemical potential 
(which may be a function of temperature), the equivalence becomes invalid due to spontaneous breaking of 
U(1)b- Our goal is now to extend the reach of the equivalence by considering deformations of the parent theory. 
Such deformations are efficaciously designed to obstruct transitions to the broken phase. 

Our motivation to consider such deformations is inspired by recent progress in a related matter, namely 
Eguchi-Kawai reduction and the volume independence of Yang-Mills at large N c , see IBhanot et al. 1982 



Eguchi and Kawai 1982 , Gonzalez-Arroyo and Okawa 1983 , which can be rephrased in the language of an 



orbifold equivalence, Kovtun et al. 2007 . Large N c volume independence relates gauge theories in different 
physical volumes provided that they remain in their center-symmetric phases. To study Yang-Mills on the 



4 At large N c , the effect of the chiral anomaly vanishes, giving one a U{2Nj) symmetry. Consequently there is one more 
Goldstone mode, which corresponds to the flavor-singlet rf meson. 

5 What is remarkable (and true of all QCD-like theories), is that in the asymptotic regime while BCS pairing is supported, there 
is no superconductor, because the BCS pairs are gauge invariant. QCD appears to be special, because BCS pairing is not gauge 
invariant. 
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lattice, it would be ideal to reduce the theory to the smallest possible volume: a single plaquette. Under this 
extreme volume reduction, center symmetry is spontaneously broken, and one can no longer stud y the physics 

that one can 



Unsal and Yaffc 



2008 



of the confined phase. To rescue volume independence, it has been argued, 
add certain double-trace terms to the original Yang-Mills action. The resulting deformed theory is the same 
as Yang-Mills in the large N c limit. The double trace terms, however, can be chosen to prevent the breaking 
of center symmetry, and enable the volume reduction to a single plaquette. This idea is being tested in lattice 



2009 



simulations, and has been met with favorable success, Bringoltz and Sharpe 

For the case at hand, we wish to stabilize SO{2N c ) gauge theory at finite density against the breaking of 
U(1)b- Baryon number breaks spontaneously as the chemical potential is increased because the baryonic pion 
energy is lowered to zero. The purpose of the deformation in this case is to lift the mass of the baryonic pion 
delaying the onset of the phase transition. If the baryonic pions were elementary fields of the theory, we could 
simply add a mass term for these modes. In terms of elementary fermion fields, however, the best we can do 
is add a four-fermion term that leads to a repulsive interaction between fermions coupled in the scalar diquark 
channel. Two particularly useful choices for the deformation are given explicitly by, Cherman et al. 2011 



Cherman and Tiburzi 2011 



N f 
a . b 



± P ab Pab 



(5) 



where S a b — ipaCjs^b has the quantum numbers of a scalar baryonic pion, and P ab = i/jaCipb has the quantum 
numbers of a pseudo-scalar baryonic meson. Because C 2 > 0, scalar diquark scattering is repulsive, and the 
mass of the Bb baryonic pions should be additively renormalized by an amount proportional to C 2 . If this is 
indeed the case, the U(1)b breaking phase transition has been offset by an amount proportional to C 2 . 

Although the deformation is physically well-motivated, the non-perturbative dynamics of the system deter- 
mine its effects. In general, the study of the deformed theory requires a lattice computation. One can use 
effective field theory techniques, however, to investigate the effects of the deformation at low-energies. When 
the fermion mass, chemical potential, and deformation are set to zero, the low-energy dynamics are described in 
terms of the emergent Goldstone modes. These modes are packaged in the coset field, £ g SU(2Nf)/ SO(2Nf). 
One then treats the chemical potential and fermion mass terms as perturbations about the chiral limit. The 
same can be done for the four-fermion operators contained in the deformation. In fact, four-fermion oper- 
ators have long been treated in chiral perturbation theory, see Bernard et al.| [1985] , A particularly simple 



way to account for such four-fermion operators is to use spurious, for a related example, see Bar et al. 2004 



Tiburzi 2005 . Using a spurion analysis, the effects of the deformation on the low-energy theory have been 



determined, ICherman and Tiburzi 2011 . For example, chiral perturbation theory for the V + deformed theory 



contains an additional term, 



Vf = v+J2 [Tr (£L (ab) ) Tr (V Z> 6 >t) + Tr (y,R^ Tr (tf R^tt 

i,6 



(0) 



where the L^ ab ^ and i?( ab ) are the spurion fields used to make the interaction symmetric under SU(2Nf). At the 
end of the day, these fields take on prescribed constant values. The parameter v + is a new low-energy constant. 
Its value is directly proportional to C 2 , and has been shown to be positive. 

Expanding the £ field out in terms of Goldstone modes, £ = 1 + i$ + . . ., the above operator gives rises to 
masses of the Bb modes. These masses are proportional to v + and so delay the onset of the U(1)b breaking 
phase transition. In fact, the vacuum structure of the V± -deformed theories has been analyzed in detail using 
the effective theory. In Fig. [2j we show the schematic phase diagrams of the deformed theories as a function of 
the baryon chemical potential squared. From the figure, one sees that the BEC phase can be avoided by tuning 
the deformation to be large enough. As a result, one can remain in the normal phase for which the orbifold 
equivalence is valid. The effective theory, hence, points to success for the deformations, however, there are a 
number of practical issues that must be resolved in order to study the deformed theories using lattice gauge 
theory techniques. We mention some of these issues in our concluding words. 



5. Summary and Outlook 



Above we have investigated the orbifold equivalence between SO(2N c ) gauge theory at finite baryon chemical 
potential and SU (N c ) gauge theory at finite baryon chemical potential. This orbifold equivalence is novel, 
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Figure 2: Schematic phase diagram of deformed SO(2N c ) gauge theory. The x-axis, where V± — 0, corresponds to 
the undeformed theory, and V± represents the strength of the deformation. Regions where baryonic-pions form a BEC 
correspond to those regions where the orbifold equivalence to large iV c QCD breaks down. Notice that one can remain 
outside the BEC phase for a large enough value of the deformation. There is a large N c exotic phase characterized by 
a BEC of baryonic pions in addition to an rj condensate. Such a phase only exists in the ^--deformed theory, and has 
been shown to be metastable. 



because the former theory does not have a sign problem at finite baryon density, while the latter does have a 
sign problem. In regimes where the orbifold equivalence is valid, one can study aspects of the sign-problematic 
theory by using lattice Monte-Carlo simulations of the sign-problem free theory. The correlation functions which 
have non-trivial images in the sign-problematic theory are identical to those in the sign-problem-free theory up 
to corrections suppressed by powers of 1 /N c . 

Using the technique of orbifold equivalence, the fermion sign problem at large N c has thus been mapped into a 
validity problem at the non-perturbative level. At relatively small values of the chemical potential, /is > |m T , 
the equivalence is known to be invalid due to spontaneous baryon number violation in the parent theory. This 
situation is analogous to that of the phase quenched approximation to QCD at finite baryon density. 6 For small 
chemical potential, it might be harmless to quench the phase of the fermion determinant, but at hb > ^rn^, 
the approximation is horrible due to the formation of a pion condensate, which has no analogue in the phase 
unquenched theory. Orbifold equivalence shows that all three theories are basically the same at large N c as long 
as fin < \m % . Beyond this value, only the phase quenched theory and SO(2N c ) theory remain equivalent. 

To extend the equivalence to larger values of the chemical potential, it is necessary to deform the parent 
theory by terms that are suppressed by powers of 1/N C . Such deformations are designed to preserve U(1)b 
symmetry as the baryon chemical potential is increased. The effects of the deformation on the dynamics of 
SO(2N c ), however, can only be determined non-perturbatively. Fortunately at low-energies, we can investigate 
the effect of the deformation using effective field theory, and the analysis indeed points to the preservation 
of U(1)b- Our effective field theory analysis assumes the power counting: \j? B ~ m 2 ~ C 2 . In this way, the 
four-fermion operators have been promoted from irrelevance to the same order as other terms, with these terms 
considered perturbatively small compared to the chiral symmetry breaking scale. To extend the equivalence to 
progressively larger values of the chemical potential, one will very likely need to increase the value of C 2 . There 
is a limit to how far C 2 can be increased in a lattice calculation, because the deformed theory formally has no 
continuum limit. It may be possible that certain aspects of the large iV c limit allow C 2 to be taken arbitrarily 
large, and this issue is under active investigation. 

A further issue of practical importance is the inclusion of four-fermion operators in the lattice theory. Because 



Recall that phase quenched QCD is identical to QCD with an isospin chemical potential. 
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fermion fluctuations are integrated out to evaluate the path integral, one must first have a Gaufiian integral for 
the fermions. The only known way to trade four-fermion terms for fermion bilinears is by using auxiliary fields. 
Auxiliary fields are integrated in using the Gaufiian integral 



1 



dye- y+Mxy . (7) 



In this toy model example, the four-fermion interaction, x 4 , is traded in for a quadratic interaction vertex, 
x 2 y, with an auxiliary field y. If one is unlucky, the introduction of auxiliary fields can reintroduce the sign 
problem. Currently it is only known how to introduce the V- deformation without a sign problem in the chiral 
limit, [Cherman et al. 2011 . An obvious avenue for future work is to devise a different implementation of the 
auxiliary fields that al ows one to simulate away from the chiral limit. 

Another important problem is to understand the mapping between baryonic observables between the SO(2N c ) 
theory and large N c QCD. Baryon interpolating operators are Z2-neutral when N c is even, suggesting that 
baryons may be in the common sector of the two theories; but, the details of applying the orbifold equivalence 
to operators involving color-space epsilon tensors have not yet been worked out. With any luck, the deformed 
SO(2N c ) gauge theory may support a new phase at finite density in a region before the equivalence breaks 
down. It would be very exciting to see this from lattice studies. 

Nonetheless, we see that the technique of orbifold equivalence has shed new light onto the very difficult 
problem of finite density QCD. We hope that further activity in this direction can help provide deeper insight 
into the fermion sign problem, and ultimately increase our understanding of the QCD phase diagram. 
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